Conventional methods fail to measure Cp(u>) of glass-forming liquids 
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The specific heat is frequency dependent in highly viscous liquids. By solving the full one- 
C ' dimensional thermo-viscoelastic problem analytically it is shown that, because of thermal expansion 

and the fact that mechanical stresses relax on the same time scale as the enthalpy relaxes, the 
plane thermal-wave method of does not measure the isobaric frequency-dependent specific heat 
Cp(u). This method rather measures a "longitudinal" frequency-dependent specific heat, a quantity 
<-h \ defined and detailed here that is in-between c p (u>) and c v (io) . This result means that no wide- 

Oh. frequency measurements of c p (uj) on liquids approaching the calorimetric glass transition exist. We 

briefly discuss consequences for experiment. 



< 

00 



C3 



T3 

o 
o 



(N 

o 



C3 



-a 



- 1—1 

x 



PACS numbers: 64.70.Pf; 65.20.- 



I. INTRODUCTION 

o . 

All liquids are viscoelastic, i.e., flow only on long time scales, whereas they are elastic on short time scales. This 
fact was first described quantitatively in 1867 by James Clerk Maxwell [![. If 77 is the (dc) shear viscosity and G^ 
the instantaneous shear modulus, the Maxwell model implies that the characteristic time r separating "long" and 
"short" times is given by 



For liquids like ambient water r is in the picosecond range; for highly viscous liquids r can be seconds, hours, or even 
months. 

When a liquid is supercooled, its viscosity increases dramatically upon continued cooling; since Goo is not nearly as 
temperature dependent, the relaxation time r also increases enormously. At some point - if crystallization is avoided 
- the liquid freezes into a glass. The glass transition takes place when the liquid is unable to fully equilibrate on 
\ the experimental time scale [2| ■ The properties of viscous liquids approaching the glass transition are far from well 
understood [1, 0, ^ B S 9 ■ A viscous liquid is characterized by its relaxation time which is the time it takes for the 
system to relax to equilibrium after an externally imposed disturbance. The disturbance may be electrical, thermal, 
or mechanical. The relaxation time depends somewhat on which quantity is measured, but as a first approximation all 
relaxation times may be regarded as identical. Since shear stresses also relax on this time scale, the liquid relaxation 
time is basically the Maxwell r of Eq. (TTJ) . Indeed, numerous experiments on viscous liquids confirm that the Maxwell 
relaxation time determines the rate of molecular motion beyond pure vibration 0, 0j B @1 • 

When probed by a linear-response experiment on time scales shorter than r, viscous liquids exhibit a solid-like 
response. As an example, consider dielectric relaxation. This is an experimental technique where accurate data 
nowadays are easily obtained over 9 or more decades [9j, [10(. Upon increase of the angular frequency lo, the dielectric 
constant e{u>) drops from its value when lot < 1 to a lower value at high frequencies (lot 3> 1) where dipolar rotations 
no longer contribute. The high-frequency dielectric constant corresponds to that observed in the glass where the 
molecules are basically frozen. Another example is the compressibility. This quantity is also frequency dependent and 
is also larger at low frequencies (lot C 1) than at high (lot ^> 1) Again, the low compressibility observed at high 
frequencies reflects glassy behavior. 

It has been known for almost a century that the specific heat is higher in the liquid phase than in the glass - indeed, 
the glass transition is often identified from specific heat measurements ;12]. It thus appears that the specific heat in 
the viscous liquid phase must be frequency dependent. This possibility was discussed briefly by Angell and Torell in 
1983 [13], but was also implicit in earlier theoretical papers [3,[lH[i3]- The frequency-dependent specific heat of a 
viscous liquid was first measured in 1985 by Birge and Nagel and by Christensen in two independent works utilizing 
different methods [IT], [l8|, [l9| . By coincidence both works measured on glycerol, but the results were not consistent 
(a discrepancy that is still unsettled, incidentally) . We here focus on the Birge- Nagel method; it covers a much wider 
frequency range than that of Ref . [lj| . 

Perhaps because of its unusual thermodynamic reference there were early doubts about the nature of the frequency- 
dependent specific heat [20I . |2ll |22| . These critiques were based on the assumption coming from generalized hydrody- 
namics that, although transport coefficients may well be frequency dependent, linear constitutive quantities do not 
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depend on frequency. This is an assumption, however, and not based on first principles. Today few researchers doubt 
that the frequency-dependent specific heat is a standard linear-response quantity; in fact, it is straightforward to derive 
a fluctuation-dissipation theorem for this quantity [23| . Thus just as for any other linear response, upon a periodic 
small time-dependent temperature variation ST(t), the enthalpy variation may be written as a convolution integral 
over the past temperature history. The frequency-dependent specific heat is the Laplace transform of the kernel of 
this convolution integral and its imaginary part determines the dissipation (i.e., free energy loss) I?! 19 L 23l|. That 
the specific heat is frequency dependent is obvious in simple energy-based master-equation models [23l . [24 1 ; in fact 
even the simple one-dimensional asymmetric double- well potential has a frequency-dependent specific heat. Recent 
works brought further insights in this continuously developing field, both as reg ards novel experimental techniques 
HE US EE EE El, El and theory including computer simulations [|1 El El, Ei • 

The method of Birge and Nagel [IE EH is based on so-called thermal effusion. A thin metal film evaporated on a 
solid glass substrate in the form of a plane slab is immersed in the liquid. The metal film acts both as heat generator 
and thermometer [IE EE EE EE EE EH ■ An applied oscillatory electric current at angular frequency uj e generates 
oscillating Joule heating at the frequency uj = 2u e due to the electric resistance of the metal film. The generated power 
"effuses" as thermal waves into the liquid and into the supporting slab. The resulting temperature oscillations depend 
on both specific heats and heat conductivities. The amplitude and phase of the temperature oscillations are detected 
by the clever so-called 3w e -technique [IE EH : Since the electric resistance of the film is temperature dependent, the 
resistance is perturbed at frequency 2u e in concert with temperature. The detected voltage across the film thus 
contains a 3w e component proportional to temperature arising from the product of varying current and resistance. 
The temperature is measured at the heat-producing film surface itself. If the thermal wavelength is short compared 
to the lateral dimensions of the film, the film is effectively an isothermal surface in space (not in time, of course). 
The heat current is orthogonal to this surface, implying that the thermal waves are planar and that boundary effects 
can be neglected. In the periodic situation all fields vary with time proportional to exp(si) where s = ±iuj depending 
on convention (4lj . By solving the heat diffusion equation Birge and Nagel arrrived at the following expression for 
the (area-specific) thermal impedance of the liquid, Z, i.e., the ratio between the complex amplitudes of temperature, 
ST(uj), and heat current density, j(u>): 



z - ^ = M)- 1/2 • (2) 

Here A is the heat conductivity and c p is the (generally frequency-dependent) specific heat per unit volume. In this 
paper we prefer to discuss the reciprocal function, the thermal admittance Y = 1 jZ which is additive when taking the 

thermal response of the substrate into account: Y = (sciXi) 1 ^ 2 + (sc 2 A 2 ) 1 ^ 2 where index 1 refers to liquid properties 
and index 2 to those of the supporting slab. Denoting the complex amplitude of the total power by P{oj), the total 
admittance Y to t — P(oj)/8T(ui) scales with the area A of the film, Y to t = YA. 

The terms "thermal impedance" and "thermal admittance" come, of course, from the close analogy between flowing 
electricity and flowing heat [i^. Temperature is analogous to potential and heat current to electric current. Thus 
Fourier's law becomes equivalent to Ohm's law, and the definition of heat capacity (specific heat) becomes equivalent 
to the definition of electric capacity. Even the dissipation of the free energy - being proportional to the real part of 
the impedance - follows from this scheme, justifying the terminology. 

According to Eq. ||2J| the bulk property measured by the Birge-Nagel experiment is the effusivity VAc EE EH • 
If the heat conductivity is frequency independent, as was generally assumed and recently confirmed [4a , 147] . the 
method yields information about the frequency-dependent specific heat. The main conclusion of the present paper, 
however, is that the experiment does not measure the isobaric specific heat. This is because the ordinary heat diffusion 
equation fails for highly viscous liquids. The precise nature of the specific heat measured can only be clarified from a 
detailed analysis. This analysis is the subject of the paper. A brief description of the problem is the following: The 
ordinary (dc) specific heat exists in isochoric, cy, and isobaric, c p , versions. Similarly, there are two ac specific heats, 
cy(w) and Cpipo). Since the pressure is usually ambient pressure, one would a priori expect that experiments measure 
c p (lo) which, indeed, is what is always reported in the literature. There is a problem, however, which implies that 
c p (uj) is not the quantity measured by the plane-plate experiment. The problem is that the coupling between the 
temperature field and the displacement field induced by thermal expansion implies that the ordinary heat diffusion 
equation fails. In other situations this can be ignored, but a subtle combination of two physical conditions prevailing 
at the glass transition renders the usual approach invalid: First, the quantity (c p — cy)/cv has a non-negligible value 
in the viscous liquid state. Secondly, upon increasing frequency the ratio of shear to bulk modulus increases to a 
non-negligible value. Physically what happens is that the high viscosity hinders the release of stresses induced by the 
thermal wave. Thermal stresses relax on the very time scale that experiments focus on. This implies that the specific 
heat is not measured at isobaric conditions; the stress tensor is not proportional to the unit tensor because the liquid 
is not in hydrostatic equilibrium. 
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For solids the difference between c p and cy is small and can usually be neglected. For less-viscous liquids internal 
stresses are quickly released and the stress tensor is proportional to the unit tensor to a good approximation. The 
problem of measuring c p (u>) reliably arises only close to the glass transition. Once the problem is recognized, the 
following questions arise: Which quantity is measured, and how should one proceed to obtain reliable c p (u>) data? 
In this paper we present the formalism needed to answer these questions by treating the one-dimensional case in 
detail. From first principles we derive the results for the adiabatic boundary condition that were simply stated in 
Ref. [48l | : these results are supplemented by results for the isothermal boundary condition. After reviewing the 
basic equations governing thermo-viscoelasticity we derive the general solution of the one-dimensional case relevant 
for understanding the results of the Birge-Nagel plane-plate method. The solution is based on the transfer matrix 
formalism, a little- known but very convenient formalism which is introduced by first discussing the simpler special case 
of no thermo-mechanical coupling. Based on the general solution the thermal admittance is worked out for several 
cases corresponding to various experimentally relevant boundary conditions (adiabatic or isothermal, mechanically 
clamped or free). Finally, the results are discussed in light of what they tell us about the prospects of obtaining 
reliable c p (u>) data by means of the plane plate method or otherwise. 



II. BASIC EQUATIONS OF THERMO-VISCOELASTICITY 



A. Thermoelasticity 

Thermoviscoelasticity describes the coupling between thermal and mechanical deviations from equilibrium. We 
shall only discuss the linear case which is well understood on the phenomenological level. To illustrate the general 
method we first discuss the simpler case where all linear- response coefficients are frequency independent. This is 
referred to as thermoelasticity, a theory that describes ordinary solids well [4£| . 

The linearity requirement implies that the system is assumed to be infinitesimally close to equilibrium. The mass 
density is denoted by p, its average by po; the temperature is T with the average temperature To. Deviations 
from equilibrium are quantified in terms of the displacement field u = u(r,i), temperature variation by <5T(r,i) — 
T(r, t) — Tq. The equations governing these variables involve the isothermal bulk modulus Kt and the shear modulus 
G (isothermal and adiabatic shear moduli are identical). Recall that these quantities are defined as follows: If V is 
volume and p pressure, the isothermal bulk modulus (inverse compressibility) is defined by 

The shear modulus G is by definition the proportionality between an off-diagonal component of the stress tensor and 
the corresponding shear deformation. 

If the heat current density is denoted by jq, the basic thermoelastic equations involve the following three constitutive 
quantities: The heat conductivity A defined via jq = — AV<5T, the isochoric specific heat per unit volume cy, and the 
isochoric "pressure-temperature coefficient" (3y defined by 



The thermoelastic "equations of motion" are nasi sua 

d 2 u / 4 



Po-q-2 = [K T + -Gj V(V • u) - GV x (V x u) - (3 V VST (5) 
d ST d 

Cy ~dT + ^ T o^(V-u) = AV 2 «\ (6) 

For completeness, we briefly sketch how these equations are derived: If di is the derivative with respect to the i'th 
spatial coordinate xt and Ui is the i'th component of the displacement vector u, the stain tensor is defined by 
€ij — (diUj + djUi)/2. Using the Einstein summation convention and denoting the stress tensor by cry , we first note 
that the relevant constitutive relation is the Duhamel-Neumann relation [50 . |51| , 

0y = K T e kk 5 i: j + 2G (eij - -tkkSij] - f3yST5ij . (7) 
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Considering a few special cases quickly convinces one that this equation is correct: If temperature is constant and 
the deformation is isotropic, only the first term contributes; thus Eq. |(7J) reproduces the definition of the isothermal 
bulk modulus because e^k — V • u is the relative volume change and the pressure is the negative diagonal element 
of the stress tensor. If temperature is constant and a pure shear deformation is considered, the second term gives 
the definition of the shear modulus. Finally, if there are no mechanical deformations, the last term reduces to the 
definition of f3y, Eq. ((4|. Now, Newton's second law for a small volume element AV is: pAVili — Fi. Here Fi is 
the i'th component of the force which is AV times the divergence of the stress tensor: Fi = AVdjdij. Equation 
((SJ) follows by substituting Eq. ((7|) into Newton's second law, because p to lowest order may be replaced by po- 
Equation (J6j> is based on the following thermodynamic relation for the volume element AV: SS = (dS/dT)yST + 
{dS/dV) T S(AV) = (AVc v /T )5T + /3 V S(AV) where the last equality follows from the definition of c v and the 
Maxwell relation (dS/dV)x = Pv- F° r the entropy per unit volume s = S/AV, since the relative volume change is 
the divergence of the displacement field, 5(AV) / AV = V • u, we get 6s = (cy /Tq)ST + /3yV • u. Entropy is conserved 
to first order: (ds/dt) + V • j s = where j s = — (A/Tq)V(5T) is the entropy current density. Thus to lowest order we 
have {ds/dt) = (c v /T )(d5T/dt) + f3yd{V ■ u)/dt = \/T V 2 (ST) which is the required Eq. ©. 

Before proceeding we note that if a p is the isobaric thermal expansion coefficient (a p = (l/V)(dV/dT) p ), the 
mathematical identity (dp / dT)y(dT / dV) p (dV / dp)r — —1 implies that 

Pv = a p K T . (8) 

In particular, if there is no thermal expansion upon heating isobarically (a p = 0), one has f3y = 0. In this case Eqs. 

and ^ decouple and reduce to the ordinary elastic equation of motion and Fourier's equation for heat conduction, 
respectively. Thus the coupling between mechanics and thermodynamics arises only when there is nonzero thermal 
expansion - which is indeed intuitively obvious. 

B. General thermoviscoelastic equations 

We now turn to the general linear thermoviscoelastic case. According to the "correspondence principle" (4^, linear 
time-dependent equations are arrived at from the thermoelastic equations ^ and ^ by replacing each time derivative 
by an s factor: 

Po s 2 u = [k t + ^G^j V(V • u) - GV x (V x u) - (3 V VST (9) 

c v sST + PvTqsV u = XV 2 5T. (10) 

In these equations the constitutive coefficients Kt, G, (3y, cy and A are generally frequency dependent, but for 
simplicity we shall not explicitly indicate this. In this way the theory is generalized from thermoelasticity to ther- 
moviscoelasticity. It is important to note that th e sy mmetry of the constitutive equations are maintained; thus the 
Maxwell relations become Onsager relations [l5|, [13, [Hj]. From here on we exclusively use the description in the 
frequency domain. 

III. FULL ANALYTIC SOLUTION OF THE ONE-DIMENSIONAL CASE 

This section presents the solution of the one-dimensional case, i.e., the case where all functions vary only in one 
spatial direction x. The transfer matrix formulation of the solution given below is convenient in practice for discussing 
various experimentally relevant boundary conditions |42j |. 

A. Case of no thermo-mechanical coupling 

As an introduction to the transfer matrix technique we first solve the case of no thermo-mechanical coupling 
(flv — 0) which, although much simpler than the general case, maintains an independent interest. If c = cy {=c p 
when j3y = 0) is the specific heat per unit volume, A the heat conductivity, and j the ^-component of the heat current 
jq, the equations for temperature and heat current are 
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ST ^ - M(x) ( Cl ) , (14) 



dST dj 

dt dx 

, dST 

In a steady-state periodic situation where all variables vary in time propotionally to exp(si) these equations become 

dj 

csST = - — 
ox 

It is convenient to discuss the problem in terms of the two variables ST and j, although one usually focuses on the 
heat-diffusion equation which results from combining the two equations: csST = Xd x 6T. The general solution of this 
equation is 8T(x) — c\ exp(fcir) + ci exp(— kx) where c\ and C2 are integration constants and 

k 2 = y ■ (13) 
Now j(x) is found from 5T(x): j(x) — — Xk[c\ exp(/ca;) — C2 exp(— kx)]. The results may be summarized as follows: 

where 

/ gkx g kx \ 

M(x) = ^_ Xkekx Xke „ kx ) . (15) 

This formulation allows one to establish the connection between temperature- and heat-current amplitudes at one 
point, say x = 0, and at the point at x: By eliminating the integration constants one arrives at 

where the "transfer matrix" T(x) is defined by 

T(x) = M(x)M(0)- 1 . (17) 
By straightforward calculation one finds that the transfer matrix is given by 

( _ ( cosh(fcx) -sinh(A;a;)/(Afc)\ , . 

[ > ~ ^-Afcsinh(fca;) cosh(fca;) ) - [ ' 

Further calculation shows that the consistency requirement T(x)T(y) = T(x + y) is fulfilled; in particular T(0) = 1 
andT(-i)=T(i)- 1 . 

From the general solution various special cases are easily derived. Suppose the sample is located between x = and 
x = L and heated periodically at x — 0. To calculate the frequency-dependent (area-specific) "thermal admittance" 
Y experienced from the heating side of the sample, Y = j(Q)/ST(0), we note that there are two experimentally 
relevant boundary conditions referring to the sample end at x = L: the adiabatic boundary condition, j(L) = 0, and 
the isothermal boundary condition, 5T(L) = 0. In the adiabatic case Eqs. (JTSJ) and (TTS|) imply that = j(L) = 
-Afcsinh(fcL)(5T(0) + cosh(fcL)j(0) which implies 



^adiab = Afctanh(fcL) 



(19) 
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On the other hand, the isothermal boundary condition at x = L is = 6T(L) = cosh(fcL)i5T(0) — sinh(fcL) J(0)/(Afe) 
which implies 

^ soth = Afccoth(fcL) . (20) 

Before proceeding to the general thermoviscoelastic case, we identify the limits of samples that are thick and thin, 
respectively, compared to the frequency-dependent thermal diffusion length: The complex k defines a (complex) 
characteristic length via In = l/k. H L ^> \ln\ one is in the "thermally thick limit" [54|, a case that applies at high 
frequencies where the thermal wave does not reach the end of the sample. In this case the two solutions become 
identical: 

*kdiab = Visoth = Afc = (sAc) 1 / 2 (\kL\ » 1). (21) 

In the "thermally thin limit" at low frequencies [54| . L <C \Id\ , the adiabatic boundary condition leads to 

^diab = Afc 2 L = scL (\kL\ <§; 1) , (22) 
corresponding to a capacitor. For the isothermal boundary condition the admittance becomes frequency independent: 



lisoth = 2 (\kL\ « 1) . (23) 

In this limit the magnitude of the heat current is solely determined by the temperature difference between the points 
x = and x = L, corresponding to the static limit of ordinary heat conduction. 



B. Full one-dimensional thermo-viscoelastic problem 



We now proceed along the above lines to determine the 4x4 transfer matrix for the full one-dimensional thermo- 
viscoelastic problem, restricting ourselves to frequencies much lower than those where sound waves become important, 
the inertia-free limit. All relevant fields depend only on the ^-coordinate and all displacement is in the x-direction. 
This does not mean that, e.g., a yy is zero - on the contrary, the stress tensor adjusts itself to be consistent with the 
clamping of all motion not in the exact x-direction. Physically, this situation corresponds to a box-shaped sample 
where the boundary conditions on the faces perpendicular to the heating plane are adiabatic and sliding. More 
realistically, the one-dimensional case to a good approximation describes the case where the dimensions of the plane 
heater are much larger than the relevant thermal wavelength. 

We still assume periodically varying fields oc exp(si). The first of the four coupled equations to be solved comes 
from Eq. ©. In the inertia- free limit the left-hand side is zero. If the ^-coordinate of the displacement vector u is 
denoted by u, since V x u = 0, this equation when integrated over x implies (where a\ is an integration constant) 

^ = idW T + ai ' (24) 

Just as with all amplitudes and coefficients here and below, the integration constants generally depend on frequency. 
The next equation is arrived at by substituting Eq. (12"4")) into Eq. (fTU|) which leads to 

A d 2 JT = sci ST + s/3 v T ai , (25) 
where the "longitudinal heat capacity" q (a name justified by the results derived below) is defined by 



In passing we note that, utilizing Eq. © and the well-known thermodynamic identities c p — cy = Tq^Kt and 
Cp/cy = Ks/Kt where K$ = —V(dp/dV)s is the adiabatic bulk modulus, the longitudinal heat capacity obeys 
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± = Ks + 4G/3 

c v K T + 4G/3 ' 1 ' 

In analogy to the standard abbreviation 7 = c p /cv (a quantity that is also generally frequency dependent) we define 

71 = — ■ (28) 
cv 

This quantity equals 7 at low frequencies where G becomes negligible compared to Kt, but generally for viscous 
liquids approaching the glass transition one has |7;| < |7|. The third equation involves the xx-component of the stress 
tensor which by Eq. ([7]) in the one-dimensional case is given by 

o xx = {K T + 4G/3) d x u ~ (3 V 5T . (29) 
The final equation is the defining equation for the time integral Sq of the heat current j: 

ls q = j . (30) 
In terms of Sq the definition of the heat-conduction coefficient A Eq. (fTTj) becomes 

Sq = --d x ST. (31) 
s 

In order to solve the four equations Eqs. (|24p . (f2l>)) . (12"9"|) and (|3"Tj) it is convenient to first rewrite them in dimen- 
sionless units. If Id is the complex number with positive real part defined by (here and henceforth q is used in the 
definition of 

1% = —, (32) 

Q S 

the (complex) dimensionless space and time variables arc defined by 

x = x/Id 

t = st. (33) 

Any function f(x,t) is made dimensionless by writing f(x,t)/fo = f(x,t). The four basic fields of the problem are 
scaled as follows: 

u = u/Id 
ST = 5T/T 

&xx = Vxx/Kt (34) 
Sq = 8q/{K T l D ). 

We finally define the following dimensionless linear-response quantities: 

c = ToCi/K T 

g = 4G/(3K T ) (35) 
a = Toctp . 



For later reference, note that in terms of these variables the following identity applies 



When Eqs. l[23]). ([25| . ([25)1 and (J2J) are rewritten as dimensionless equations, the result is 



u = ST + a% 

1+5 

ST" = ST + ~ ai 
c 

5-3,3, = (1 + g)u - a ST 
Sq = —cST'. 

The integration of these equations is straightforward. Substituting Eq. ([57)) into Eq. ([55)) leads to 

&xx = (1 + g) a\ ■ 

The integration of Eq. (|38p introduces two new (frequency-dependent) integration constants: 



ST = e x a 3 + e * a 4 - — a x 

c 

This is substituted into Eq. (J37j, the integration of which results in 



u = 1 



ae ae 
x oi + a 2 + a 3 — - i _ a 4 



2(1 + 5),/ ^ ' * ' 1+5 " 1 + 
Finally, when Eq. (jHZ)) is substituted into Eq. ([¥0)1 the result is 



Sq = —ce x a 3 + ce x a 4 . 
The complete solution is summarized as follows 



ST 
u 

\SqJ 



M(f) 



where 



M(f) 



/ 1+3 \ 

-S r r e" £ 



x 1 



ae ae 



V -ce x ce~ x J 



The transfer matrix T(f ) connecting the fields at x — with those at x, 



( o xx (x)\ 
ST(x) 
u(x) 

V Sq(x) J 



T(£) 



/?x«(0)\ 

<5T(0) 
u(0) 
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is given by Eq. (fl7|) . The explicit calculation of T(f ) results in 



T(f) 





cosh(i) 

1 



/ 1 

a(cosh(x) — 1) 

WW 

a (sinh(x) — x)-\-c(l-\-g)x a sirthfaf) 

W+FF i+s A WW 

_asn^x)_ — csinh(a;) cosh(J) ) 



° \ 

sinh(a;) 

c 

a(l — cosh(x)) 



(48) 



As in the case of simple heat conduction, a direct calculation shows that the consistency requirement T(x)T(y) = 
T(x + y) is fulfilled, implying that T(0) = 1 and T(-x) = T(f) -1 . 

IV. THE THERMAL RESPONSE FOR A NUMBER OF BOUNDARY CONDITIONS 

Consider now a situation where at the surface x — one imposes a heat current density with amplitude j(0) and 
measures the corresponding temperature amplitude 6T(0). The liquid is assumed to be fixed at the x — surface, 
i.e., u(0) = 0. The experimentally relevant quantity is the thermal (area-specific) admittance, Y — j(Q)/ST(0). If one 
defines the dimensionless thermal compliance by 



J Sq(0) 



in real units one has 



Y 



ST(0) 



K T l D s - 
To J - 



We proceed to calculate Y for a number of cases by means of the inverse transfer matrix, 



(49) 



(50) 



P = T" X (L) = T(-L) . 
For reference we give the explicit matrix P found by substituting x = — L into Eq. (|48|): 



(51) 



/ 1 

q(cosh(£)-l) 
W+t) 

q 2 (- sinh(L)+L)-c(l+g)L 
c(l+g)2 
a sinh(L) 
1+9 





cosh(Z) 



a, sinh(L) 

i+l - 



o 

sinh(L) 

c 

q(l-cosh(L)) 



WW) 

csinh(L) cosh(L) J 



(52) 



For each case we shall explicitly give the low frequency (thermally thin) and high frequency (thermally thick) limits. 

Case A. Free mechanical termination, a xx {L) = 0. 
Subcase 1. Adiabatic termination. 

This case corresponds to Sq(L) — and thus Eq. (|47|) implies 



(53) 







( °A 


Sf(0) 


= P 


ST(L) 







u(L) 


V Sq(0) J 




\ o j 



From this we get 



J = 



m = P i2 6f(L) + P i3 u(L) _ P42 _ L 
6T(0) P 22 6T(L) + P 23 u{L) P22 



(54) 
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Thus the thermal admittance is given by 



Y = cilps tanh(L/?£)) = y / ACista,nh(L/lo) ■ (55) 

In the thermally thin limit, L <C \Id\, this reduces to Y = \JXcisLjlu — ciLs. In the thermally thick limit, L> \^d\, 
one finds Y = \J Xcis. 

Subcase 2. Isothermal termination. 
This case corresponds to 5T(L) — and thus Eq. (|T7)l implies 







< \ 


Sf(0) 


= p 










u(L) 


V Sq(0) J 




\Sq(L)J 



(56) 



From this we get 

J= m = WM| = ^ = 5cotn(£) . (57) 

ST{0) P 23 u{L) + P 2i 5q{L) P2A 
Thus the thermal admittance is given by 

Y = c 4 l D s coth{L/l D ) = y/Xc^s coth(L/l D ) . (58) 

In the thermally thin limit, L <C \Id\, this reduces to Y = V Xcislo/L = X/L. In the thermally thick limit, 
L S> \Id\, one finds Y = ^/Xcjs. 

Case B. Clamped mechanical termination, u(L) = 0. 
Subcase 1. Adiabatic termination. 

This case corresponds to Sq(L) = and thus Eq. (jl?)) implies 



fa xx (Q)\ 




/ a xx (L)\ 


Sf(Q) 


= P 


5f(L) 










V Sq(0) J 




V o / 



(59) 



From the second and fourth equation of Eq. ([59")) we get 



j = 5q{fi)_ = P 41 +P 42 Sf(L)/a xx (L) 
Sf{0) P2i+P225T{l)/5 xx {L)' 



and from the third we get = P-nd xx {L) + P^2&T{L) or 



Sf(L)/a xx (L) = (61) 

^32 



Substituting this into Eq. ((60)) leads to 



A calculation of this leads to 



P ^ - P «*" . (62) 
P21P32 - P22P31 



J = = . (63) 

coth(L) + ( 7I - 1)/ L 
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Thus in real units I5E 



Y = 



cilos 



(64) 



coth(L// D ) + {l D /L){ci/c v - 1) ' 

In the thermally thin limit, L <C \Id\, this reduces to Y = cyLs. In the thermally thick limit, L S> |Zd|, one finds 
Y = ^/Xc7s. 

Subcase 2. Isothermal termination. 
This case corresponds to 5T{L) — and thus Eq. (|37[) implies 







/ a xx (L)\ 


Sf(0) 


= p 













V Sq(0) J 




I Sq(L) J 



(65) 



From the second and fourth equation of Eq. (|65|) we get 



J = 



Sq(0) = PAi+Pa5q{L)/a xx {L) 
6f(0) P 2 i+P24Sq(L)/a xx (L) 



(66) 



and from the third we get 



Sq(L)/a xx (L) 



P31 

>34 



(67) 



Substituting this into Eq. (f6"6")) leads to 



J = 



PaPzi — -P44-P31 
P21P7H — P24P31 



(68) 



or 



J = 



(7; - 1) sinh(Z) + icosh(Z) 
2(7; - l)(cosh(L) - 1) + L sinh(L) 



(69) 



Thus in real units 



Y = ciIds 



(70) 



{ci/c v - 1) sinh(L/^ D ) + (L/l D ) cosh(L/l D ) 
2( Cl /c v - l)(cosh(L/Z D ) - 1) + (L/l D )smh{L/l D ) ' 

X/L. In the thermally thick limit, L \Id\, one finds 



In the thermally thin limit, L <C \Id\, this reduces to Y 
Y = VAqs. 

In the thermally thin limit (|/u(w)| S> L) the above results are summarized in table U for the four combinations 
of thermal and mechanical boundary conditions at the surface x = L. It was to be expected that the isothermal 
boundary condition gives the heat conductivity in both mechanical cases, because heat conductivity does not depend 
on whether constant volume or constant stress boundary conditions apply. For the adiabatic boundary condition it 
is also expected that the admittance is related to the constant volume specific heat cy. That the admittance for free 
mechanical boundary condition is related to q and not c p is not trivial; it relates to the fact that the displacement 
associated with the temperature oscillation is longitudinal. In the thermally thick limit, |Zd(u;)| <C L we se that the 
thermal admittance becomes Y — (sqX) 1 ^ 2 and thus is related to the longitudinal specific heat for all combinations 
of boundary conditions considered (table H3J . 
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Vthin 


adiabatic isothermal 


free 
clamped 


CiLs \/L 
cvLs \/L 



TABLE I: The (area-specific) thermal admittance Y at the boundary x = for different boundary conditions at x — L in the 
thermally thin limit, |iz)(w)| 3> L. 





adiabatic isothermal 


free 
clamped 


( S qA) 1/2 ( SCi A) 1/2 
( SCi A) 1/2 (sctX) 1 ' 2 



TABLE II: The (area-specific) thermal admittance Y at the boundary x = for different boundary conditions at x = L in the 
thermally thick limit (|Z_d(o;)| L. 



V. DISCUSSION 



How does the result Y — (sqA) 1 / 2 apply to real effusivity measurements? Are the adiabatic and sliding boundary 
conditions on faces perpendicular to the heating plane a good approximation? In actual experiments the heating plate 
width W fulfils |£d(o;)| <§C W to assure one-dimensional heat flow orthogonal to the heating plane. The associated 
displacement field has a range of |Z.d(u;)|, and if the solid slab supporting the heating film is much thicker than this 
length, it will clamp motion along the plane and only allow displacement orthogonal to the plate. In this case the 
above solution is a good approximation. 

One might expect that clamping problems could be circumvented by using an unsupported film free to expand 
laterally, but that is not correct. Suppose that the film were unsupported and consequently did not stress the liquid. 
First, we note that there would be a technical problem with the 3w detection technique since thermal expansion 
straining the film would give a spurious resistance change on top of the thermally induced resistance change. That 
problem may be solved, however, and is not the main issue here. Even in this case the thermomechanical coupling, 
however, would introduce anisotropic stresses: Suppose the strain e yy were homogeneous on any plane at right angle to 
the x-axis. Then the displacement in the y-direction at the boundary of the plate would be proportional to the lateral 
dimension W , but at a distance of the order | Z x> (a;) | in the x-direction it would be zero. This would produce enormous 
shear strains in the liquid (of order e yy ry~7~ji ), the relaxation of which is controlled by the (high) shear modulus. The 

condition |Zd(o>)| W needed for one-dimensional heat flow thus simultaneously gives rise to a situation of partial 
and anisotropic clamping of the liquid. To determine precisely which effective frequency-dependent specific heat is 
actually measured in this (hypothetical) case of a film free to expand laterally requires a detailed calculation; it will 
not be simply c p . 

The "adiabatic" method of Christensen's 1985 paper [l^| is quite different from the Birge-Nagel method, but suffers 
from the same basic problem that the liquid upon thermally expanding must flow to release the stresses, and that this 
flow is partly inhibited by the very large viscosity. Thus as for most other calorimeters, the adiabatic method also 
has uncontrolled mechanical stresses and a stress tensor that is not diagonal. This means that c p (uj) is not measured. 
As mentioned, to determine precisely which quantity a given method measures requires a detailed analysis taking 
into account the properties of the calorimeter materials. It is quite possible that the discrepancy between the 1985 
Birge-Nagel and Christensen measurements on glycerol derives from the two methods de facto not measuring the same 
quantity. 

Studying the frequency-dependent specific heat via effusion can be done in other geometries than the planar. Birge 
and Nagel (l8j also made experiments in an axially symmetric geometry. This was done in order to separate out 
any possible frequency dependence of the heat conductivity. These experiments were not presented as a function 
of frequency. Rather, temperature was scanned for fixed frequencies which makes the analysis somewhat indirect 
(the thermomechanical coupling was not considered in this early experiment, either). Christensen and Olsen [56j 
performed an experiment in spherical symmetry and did take into account the thermomechanical coupling. A small 
drop of 1,2,4-butanetriol of radius r 2 — 0.35mm was placed around a temperature-dependent resistor (thermistor) of 
radius n = 0.15 mm. By using the 3<x>-technique the thermal compliance J tot — Y to t/s was found. In the thermally 
thin limit one has [5t| 



c v V 



(r 2 /n) 3 K s + 4/3G 
{r 2 /n) 3 K T + 4/3G 



(71) 
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where V is the liquid volume. This means that the longitudinal specific heat is measured in the limit r<ijr\ — > 1 
(thin liquid shell) and that the isobaric specific heat is measured in the limit r 2 /r 1 — > oo (thick liquid shell). These 
conditions could only be met at very low frequencies, however ( < 40mHz). 

How large is the difference between c p and c{l From the two identities for 7 and 7/ 



c v 



Ks_ 

K T 



(72) 



and 



Q M s 

li = — = -JT ( 73 ) 
cy Mt 



the following useful result follows 



c P ~ ci = 4 G c p - c v 

c p 3 Mt Cp 1 ' 

For glycerol the relative change in measured specific heat (45% at the glass transition) gives a rough upper estimate 
of (c p — Cv)/c p . The infinite- frequency shear and bulk moduli are 4 GPa and 11 GPa respectively. Thus the factor 
AG/3Mt is 33% as an upper estimate. In the case of glycerol the deviation of c/ from c p becomes 15% or one third 
of the relaxation strength. For liquids with larger expansion coefficient, the effect is larger. Note that this is a 
frequency-dependent correction and that q does not just scale with c p - it has a different frequency dependence. 
Equation (|T4[) expresses the influence of the thermo-mcchanical coupling on ci in a nutshell: At low frequencies (or 
high temperatures) G/M T is small, and q becomes c p . At high frequencies (or low temperatures) (c p — c v )/c p is 
negligible, and cj again becomes c p . This is why the problem treated in this paper is not encountered for specific- 
heat measurements on solids or on ordinary liquids. The effect is only significant at the glass transition where upon 
increasing frequency the first factor [G/Mt] increases while the second [(c p — cy)/c p ] decreases. 

Finally, we wish to emphasize that the effect discussed in this paper has nothing to do with the liquid falling 
out of equilibrium as the glass transition is approached; the effect is present in the equilibrium liquid well above the 
calorimetric glass transition temperature. Thus, even if effusivity measurements are performed in, say, the kHz regime, 
the effect is present and must be taken into account whenever the alpha relaxation time is of order milliseconds or 
longer. 



VI. CONCLUSIONS 



Fundamental physical facts make it difficult to measure the frequency-dependent isobaric specific heat for liquids 
approaching the glass transition because the temperature perturbation of a thermal experiment induces stresses. The 
high shear modulus slows down the relaxation of the deviatoric part of the stress tensor, giving rise to non-isobaric 
conditions. In planar geometry this is a problem whether experiments are done in the thermally thin or the thermally 
thick limits. 

In this paper the one-dimensional thermoviscoelastic problem was solved in terms of a transfer matrix, making it 
easy find the thermal response for any types of boundary conditions. It was found that in the thermally thick case 
the isobaric specific heat should be replaced by the longitudinal specific heat. This result applies to real effusivity 
measurement to the extent they may be considered one-dimensional. 

In conclusion, wide-frequency specific-heat spectroscopy based on thermal-wave effusion determines the longitudinal 
specific heat. The isobaric and isochoric specific heat can be calculated subsequently if the frequency-dependent 
mechanical moduli are known. - Although we here only studied the thermoviscoelastic problem in the frequency 
domain, these results also give rise to serious questions regarding the interpretation of time-domain enthalpy- relaxation 
experiments. 



14 



Acknowledgments 



This work was supported by the Danish National Research Foundation's (DNRF) centre for viscous liquid dynamics 
"Glass and Time." 



J. C. Maxwell, Phil. Trans. R. Soc. London 157, 49 (1867). 
W. Kauzmann, Chem. Rev. 43, 219 (1948). 

G. Harrison, The Dynamic Properties of Supercooled Liquids (Academic, New York, 1976). 

S. Brawer, Relaxation in Viscous Liquids and Glasses (American Ceramic Society, Columbus, Ohio, 1985). 
C. A. Angell, J. Non-Cryst. Solids 131, 13 (1991). 

C. A. Angell, K. L. Ngai, G. B. McKenna, P. F. McMillan, and S. W. Martin, J. Appl. Phys. 88, 3113 (2000). 
S. P. Das, Rev. Mod. Phys. 76, 785 (2004). 
J. C. Dyre, Rev. Mod. Phys. 78, 953 (2006). 

P. Lunkenheimer, U. Schneider, R. Brand, and A. Loidl, Contemp. Phys. 41, 15 (2000). 

F. Kremer and A. Schonhals (Eds.) Broadband Dielectric Spectroscopy (Springer, Berlin, 2002). 

T. Christensen and N. B. Olsen, Phys. Rev. B 49, 15396 (1994). 

I. M. Hodge, J. Non-Cryst. Solids 169, 211 (1994). 

C. A. Angell and L. M. Torell, J. Chem. Phys. 78, 937 (1983). 

R.-J. Roe, J. Appl. Phys 48, 4085 (1977). 

C. T. Moynihan and P. K. Gupta, J. Non-Cryst Solids 29, 143 (1978). 

C. T. Moynihan and A. V. Lesikar, Ann. N.Y. Acad. Sci. 371, 151 (1981). 
N. O. Birge and S. R. Nagel, Phys. Rev. Lett. 54, 2674 (1985). 

N. O. Birge, Phys. Rev. B 34, 1631 (1986). 

T. Christensen, J. Physique Colloq. 46, C8-635 (1985). 

D. W. Oxtoby, J. Chem. Phys. 85, 1549 (1986). 
R. Zwanzig, J. Chem. Phys. 88, 5831 (1988). 

J. Jackie, Physica A 162, 377 (1990). 

J. K. Nielsen and J. C. Dyre, Phys. Rev. B 54, 15754 (1996). 
J. Bisquert, Am. J. Phys. 73, 735 (2005). 

D. M. Paolucci and K. A. Nelson, J. Chem. Phys. 112, 6725 (2000). 

H. Huth, M. Beiner, and E. Donth, Phys. Rev. B 61, 15092 (2000). 
Y. H. Jeong YH, Thermochim. Acta 377, 1 (2001). 

L. Carpentier, O. Bustin, and M. Descamps, J. Phys. D - Appl. Phys. 35, 402 (2002). 

H. Honda, H. Ogura, and A. Chiba, Measurement Science and Technol. 13, 1437 (2002). 

A. Maesono, Y. Takasaki, Y. Maeda, R. P. Tye, J. Morikawa, and T. Hashimoto, High Temperatures - High Pressures 34, 
127 (2002). 

U. Harbola and S. P. Das, Phys. Rev. E 64, Art. No. 046122 (2001). 

P. Scheidler, W. Kob, A. Latz, J. Horbach, and K. Binder, Phys. Rev. B 63, Art. No. 104204 (2001). 

C. C. Yu and H. M. Carruzzo, Phys. Rev. E69, Art. No. 051201 (2004). 

D. Chakrabarti and B. Bagchi, J. Chem. Phys. 122, Art. No. 014501 (2005). 
N. O. Birge and S. R. Nagel, Rev. Sci. Instrum. 58, 1464 (1987). 

U. G. Jonsson and O. Andersson, Measurement Sci. and Technol. 9, 1873 (1998). 

E. Gmelin, Thermochim. Acta 304, 1 (1997). 
N. Menon, J. Chem. Phys. 105, 5246 (1996). 

I. K. Moon, Y. H. Jeong, and S. I. Kwun, Rev. Sci. Instrum. 67, 29 (1996). 
N. O. Birge, P. K. Dixon, and N. Menon, Thermochim. Acta 304, 51 (1997). 

It is a matter of convention whether one uses exp(iwi) or exp(— iuit) for the complex time dependence. The form corre- 
sponding to s — iui (common, e.g., in electrical engineering literature) is convenient because it saves a lot of minus's. The 
form corresponding to s — —iu), on the other hand, is preferable when there are plane waves because when k > a wave 
of the form exp[i(fc:r — u)t)] moves in the positive x-direction. 

H. S. Carslaw and J. C. Jaeger, Conduction of heat in solids (Clarendon Press, Oxford, 1959). 
A. A. Minakov, S. A. Adamovsky, and C. Schick, Thermochim. Acta 403, 89 (2003). 
E. H. Bentefour, C. Glorieux, M. Chirtoc, and J. Thoen, J. Chem. Phys. 120, 3726 (2004). 
Y. Kraftmakher, Phys. Rep. 356, 1 (2002). 

A. A. Minakov, S. A. Adamovsky, and C Schick, Thermochim. Acta 377, 173 (2001). 

E. H. Bentefour, C. Glorieux, M. Chirtoc, and J. Thoen, J. Appl. Phys. 93, 9610 (2003). 

T. Christensen and N. B. Olsen, Prog. Theor. Phys., Suppl. 126, 273 (1997). 

R. M. Christensen, Theory of viscoelasticity, 2nd edition (Academic Press, New York, 1982). 

L. D. Landau and E. M. Lifshitz, Theory of Elasticity 3rd Ed. (Pergamon, London, 1986). 

W. Nowacki, Thermoviscoelasticity, 2nd edition (Pergamon Press, 1986). 



15 



[52] J. Meixner and H. G. Reik, in Principen der Thermodynamik und Statistik (Springer, Berlin, 1959), Ed. S. Fliigge, Handbuch 

der Physik, Vol. 3, pp. 413. 
[53] J. I. Berg and A. R. Cooper, J. Chem. Phys. 68, 4481 (1986). 
[54] J. Caerels, C. Glorieux, and J. Thoen, Rev. Sci. Instrum. 69, 2452 (1998). 

[55] This equation appeared previously in Ref. [48fl. but with a typo: The final term missed the "-1". 
[56] T. Christensen and N. B. Olsen, J. Non-Cryst. Solids 235, 273 (1998). 



